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Abstract 

The new CRISP filter at the Swedish 1-m Solar Telescope provides opportuni- 
ties for observing the solar atmosphere with unprecedented spatial resolution 
and cadence. In order to benefit from the high quality of observational data 
from this instrument, we have developed methods for calibrating and restoring 
polarized Stokes images, obtained at optical and near infrared wavelengths, 
taking into account field-of-view variations of the filter properties. 

In order to facilitate velocity measurements, a time series from a 3D hy- 
drodynamical granulation simulation is used to compute quiet Sun spectral 
line profiles at different heliocentric angles. The synthetic line profiles, with 
their convective blueshifts, can be used as absolute references for line-of-sight 
velocities. 

Observations of the Ca II 8542 A line are used to study magnetic fields in 
chromospheric fibrils. The line wings show the granulation pattern at mid- 
photospheric heights whereas the overlying chromosphere is seen in the core 
of the line. Using full Stokes data, we have attempted to observationally verify 
the alignment of chromospheric fibrils with the magnetic field. Our results 
suggest that in most cases fibrils are aligned along the magnetic field direction, 
but we also find examples where this is not the case. 

Detailed interpretation of Stokes data from spectral lines formed in the 
chromospheric data can be made using non-LTE inversion codes. For the first 
time, we use a realistic 3D MHD chromospheric simulation of the quiet Sun 
to assess how well NLTE inversions recover physical quantities from spec- 
tropolarimetric observations of Call 8542 A. We demonstrate that inversions 
provide realistic estimates of depth- averaged quantities in the chromosphere, 
although high spectral resolution and high sensitivity are needed to measure 
quiet Sun chromospheric magnetic fields. 
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1 Introduction 



The solar atmosphere constitutes a remarkably complex astrophysical labo- 
ratory that continuously performs experiments for us to observe. One reason 
for trying to understand this complicated environment is to establish with high 
precision simple but fundamental properties of the Sun. An important example 
is its chemical composition which can be put in context with our understand- 
ing of the astrophysical processes in the interior of the Sun and other stars, 
in the Galaxy, and in the early universe. To accomplish this, we need to take 
into account the dynamics of the solar photosphere as well as the physical 
processes involved in the formation of the spectral lines from which chemical 
abundance ratios are determined. This represents a major challenge and our 
confidence in the results rely heavily on the accuracy of measurements made 
with modern solar telescopes. 

A second reason for studying the solar atmosphere, and one that is even 
more relevant in the present thesis, is to understand the mechanisms that gen- 
erate the observed fine structure, dynamics, and magnetic field and to carry 
over that understanding to other astrophysical plasmas, including other stellar 
atmospheres. Dramatic progress in this field has in recent years been made in 
part by improved theoretical simulations and in part by new solar telescopes, 
equipped with powerful instrumentation, on the ground and in space. Both 
simulations and observations clearly demonstrate that much dynamics occur at 
very small spatial scales. Obtaining accurate quantitative information that will 
allow us to confirm or refute new models requires pushing existing telescopes 
to their diffraction limit and designing future telescopes with improved spatial 
resolution, better signal-to-noise and equipped with multiple instruments to 
simultaneously diagnose different layers of the solar atmosphere. In addition, 
sophisticated post-processing methods are needed to enhance the fidelity of 
the observed data, by developing accurate methods for calibrations and for 
removing contamination from stray light due to limitations set by the Earth's 
atmosphere and optical imperfections in the telescope or its instrumentation. 
This thesis deals with the challenges of accurate measurements of quantities 
relevant to small-scale dynamics, based on observations with a major solar 
telescope: the Swedish 1-m Solar Telescope (SST) on La Palma. 

Our observational data are from two distinct, but physically connected, lay- 
ers of the solar atmosphere: the photosphere and the chromosphere. The dy- 
namics and morphology of these two atmospheric layers are completely dif- 
ferent. To a large extent, these differences can be attributed to magnetic fields: 



2 



Introduction 



whereas the gas pressure falls of exponentially with height and is roughly 10 5 
times smaller in the chromosphere than in the photosphere, the magnetic field 
strength falls off much slower. The relative importance of forces associated 
with gas pressure and magnetic field can be estimated from the ratio of gas 
pressure (P g ) to magnetic pressure (Pg), the plasma-beta parameter defined 
by: 




In the photosphere, j3 is much larger than unity everywhere, except in sunspots 
and other (mostly small-scale) concentrations of magnetic field. The photo- 
sphere is dominated by a convective energy flux, peaking just below the vis- 
ible surface. Key questions today are to understand how this energy flux is 
maintained within magnetic structures and a major challenge is to identify 
and measure the velocity signatures of any convective flows present. These 
signatures are both weak and and small-scale and their identification relies on 
whether we can establish an absolute reference for measured Doppler veloci- 
ties on the Sun. The first part of the present thesis deals with this problem. 

The second part deals with observations of the chromosphere. Here, 
magnetic fields are much weaker than in the strongest magnetic structures 
seen in the photosphere but the gas pressure is even lower. In the upper 
chromosphere, j3 < 1 thus the magnetic forces are dominant. This work 
aims at investigating the potential for diagnosing the weak chromospheric 
magnetic fields using Stokes polarimetry and sophisticated inversion 
techniques. 

To set the context, we proceed with an overview of some photospheric and 
chromospheric topics. 



1.1 The Photosphere 

The visible surface of the Sun corresponds to the photosphere, a thin layer 
of about 500 km located on top of the convection zone, where the plasma 
changes from completely opaque to transparent (Stix 2002). 



1.1.1 Granulation 

Outside active regions with strong magnetic fields, the photosphere is domi- 
nated by a dynamic pattern of bright granules surrounded by dark intergran- 
ular lanes (see Fig. 1.1). The flow in a granule resembles that of a fountain, 
where the hot plasma moves upwards inside the granules and then flows out 
towards the edge, where the cooler plasma merges with material from neigh- 
bouring granules. Gravity and pressure increase at the edge of the granules, 
accelerating the fluid downwards (Stein & Nordlund 1998). Regular granules 
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have a typical size of the order of 1 Mm, and a characteristic life time of 6 
minutes (Bahng & Schwarzschild 1961). 

Convective motions leave strong fingerprints on any line formed in the 
photosphere. An important diagnostic is the C-shaped bisector obtained from 
spatially-averaged line profiles. This effect is produced by the statistical av- 
erage of bright blueshifted profiles from granules with dark redshifted pro- 
files originating in the intergranular lanes (Dravins et al. 1981). This intensity 
weighted average is blueshifted as upflows are more heavily weighted by be- 
ing brighter and covering a larger area than the narrower intergranular lanes. 
This shift is commonly known as the convective blueshift. 

The thermodynamic history of fluid elements rising through the photo- 
sphere is described in detail by Cheung et al. (2007). The temperature decrease 
of a fluid element that moves up in the convection zone, is mostly produced by 
adiabatic expansion. As the fluid reaches the photosphere, the opacity drops 
and the fluid rapidly loses entropy by radiative cooling. At this point, the fluid 
cell is overshooting into the stably stratified photosphere, and it still interacts 
with the surroundings because it is not completely transparent to radiation. 
Fig. 1.2 shows the trajectory described by tracer fluid elements that enter the 
photosphere. The color coding indicates the sign of 2 ra( j, the heat exchange 
with the surroundings, in dark for radiative loses (2 ra( j < 0) and in light grey 
where the fluid elements are being heated (<2 ra( j > 0). Interestingly, there is no 
direct correlation between heat exchange and the temperature variation of the 
fluid elements, which is determined by a balance between adiabatic expansion 
and (non-adiabatic) heat exchange with the surroundings. 

In recent years, efforts to obtain refined estimates of solar abundances have 
given rise to some controversies (see e.g. Asplund et al. 2000). This debate 
stimulated improvements of treatment of energy transfer by radiation in 3D 
hydrodynamic simulations, in particular for the mid and high photosphere, 
where spectral lines are formed. These controversies have also stimulated 
the development of improved non-LTE calculations of spectral lines used 
for abundance calculations (Shchukina & Trujillo Bueno 2001). As a result, 
3D MHD simulations of more complex solar atmosphere dynamics involving 
magnetic fields can now also be made with improved energy transfer than just 
a few years ago. 



1.1.2 Sunspots 

The structure and dynamics of sunspots remain some of the most controversial 
topics in the solar photosphere. Sunspots constitute strong magnetic field con- 
centrations that appear in the atmosphere and present typical sizes of 12000 
km. As a first approximation, we can assume that the magnetic field acts on 
the gas in the form of a magnetic pressure, 

B 2 
2jUo 
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Figure 1.1: The photosphere imaged in the G-band at 430 nm with the Swedish 1-m 
Solar telescope. Image courtesy of M. van Noort & L. Rouppe van der Voort (ITA- 
UiO). 
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Figure 1.2: Trajectories of tracer fluid elements above a granule. The grayscale image 
show the vertical velocities at at Z = when the tracers were released. The colors in 
the trajectories correspond to Q m( \ < (dark-grey) and 2 ra( j > (light-grey). From 
Cheung et al. (2007) (reproduced with permission of the authors). 



Horizontal force balance then dictates that the sum of the gas pressure and 
magnetic pressure must be the same inside the sunspot as outside. This implies 
that 

B? L _ 

where P s is the gas pressure inside the spot and P qs is the gas pressure in the 
surrounding (non-magnetic) quiet sun. An immediate consequence of this is 
that the gas pressure inside the spot must be lower than outside and that there- 
fore the atmosphere is more transparent inside sunspots, allowing observers 
to see deeper layers of the atmosphere than in quiet Sun observations. This 
is generally known as the Wilson depression (~ 500 km), discovered by Wil- 
son & Maskelyne (1774). At the same time the H opacity decreases with 
temperature, and sunspots are cooler than their environment so the opacity 
decreases even more. Sunspots are darker than their surrounding granulation 
because convection is suppressed by the strong magnetic field, thus sunspots 
are cooler as an effect of inefficient heat transfer. 

During the past decade, the scientific debate has focused on the dynam- 
ics and structure of the penumbra and explaining why the penumbra is as 
bright as it is, about 75% of the surrounding quiet sun. Recent advances in in- 
strumentation have unveiled very fine structure in sunspots, especially in the 
penumbra. Fig. 1.3 illustrates in great detail the fine structure of the penum- 
bra of a sunspot. The blown-up section shows dark cored penumbral filaments 
(Scharmer et al. 2002) and the inner umbra. The following theoretical frame- 
works (Scharmer 2008) have become popular because they can partially re- 
produce the features observed in sunspots, in spite of representing different 
physical concepts. 
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Figure 1.3: Surroundings of a sunspot (top) and close-up view of penumbral filaments 
(bottom). The units of the axis are given in Mm. Continuum observations at 396 nm, 
by Vasco Henriques (ISP-KVA). 



1. The uncombed penumbra model (Solanki & Montavon 1993) postulates 
the existence of discrete flux tubes with homogeneous magnetic field inside 
that discontinuously changes at the boundary. Those nearly horizontal flux 
tubes are embedded in a more vertical magnetic field. This model was able 
to reproduce strongly asymmetric Stokes V profiles observed on the limb 
side of the penumbra in observations off disk center (Sanchez Almeida & 
Lites 1992). 

2. Siphon flow models (Meyer & Schmidt 1968; Montesinos & Thomas 
1997) are based on the idea that a difference in field strength between the 
two footpoints of a flux tube leads to a difference in gas pressure, driving a 
plasma flow in the direction of the footpoint with the highest field strength 
(thus lower gas pressure). Evershed flows are assumed to be steady flows 
between two footpoints with different magnetic field strengths (e.g., 
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Westendorp Plaza et al. 1997). However, these models do not explain the 
mechanism involved in such field strength difference at the footpoints. 

3. Convection and downward pumping of magnetic flux are ingredients 
added to the siphon models. As siphon models present a stationary solu- 
tion, time variations are explained by external mechanisms to the penum- 
bra. In this context, moving penumbral grains are assumed to be produced 
by a moving convective pattern in the bright side of the penumbra. Thomas 
et al. (2002) and Weiss et al. (2004) attribute the submergence of the flux 
tubes at the boundary of the penumbra to downward pumping produced 
by convective motions. Furthermore, they attribute the whole filamentary 
structure of the penumbra to the same downward pumping mechanism, ex- 
plaining the structure inside the sunspot based on mechanisms that take 
place outside. 

4. The convective gap model proposed by Scharmer & Spruit (2006). In this 
framework, the proposed mechanism that generates penumbral filaments 
is convection, in radially aligned, nearly field free gaps. The strong field 
gradients that are necessary to reproduce the asymmetric Stokes V profiles 
reported by Sanchez Almeida & Lites (1992) are assumed to be produced 
by the topology of nearly field-free gaps combined with line-of-sight gra- 
dients in the flow velocity. The Evershed flow is explained as representing 
the horizontal component of this convection. 

1.1.3 Oscillations in the solar atmosphere 

Solar oscillations were discovered by Leighton et al. (1962) with a simple 
observational technique: two simultaneous images were recorded in the blue 
and the red wings of a spectral line, respectively, and then subtracted. The 
resulting image contained intensity variations produced by the Doppler shift 
of the line. Kahn (1961) proposed that oscillations are sound waves trapped 
in the solar atmosphere. Towards the solar interior, the temperature and speed 
of sound increase with the variation of the refractive index, caused by the 
increased density. The wave is refracted until it starts to propagate upwards. 
The same process occurs above the photosphere where the waves are refracted 
back into the inner atmosphere. Observations contain an overlap of hundreds 
of modes of oscillation that effectively reach different depths. The oscillations 
in the photosphere typically have a 5-minute period and an amplitude around 
1 kms- 1 (Stix 2002). 



1.2 The chromosphere 

The chromosphere represents many challenges for solar physicists. Despite 
important discoveries during the past decade, it still remains unexplored to a 
large extent. From an observational point of view, only a few spectral lines 
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Figure 1.4: Spicules at the limb seen at H I 6563 A. Intensity gradients towards the 
limb have been filtered to enhance the off-limb part. Image taken at the Swedish 1-m 
Solar Telescope, courtesy of Luc Rouppe van der Voort (ITA-UiO). 



are sensitive to the chromospheric height range and those are usually hard to 
model and understand, like H 1 6563 A, Ca n K & H (3934 and 3968 A respec- 
tively), the Call infrared triplet (8498, 8542 & 8662 A), Nal Dj (5896 A) , 
He I 10830 A. 

During the past 15 years combined efforts from observational and com- 
putational approaches have lead to a better understanding of chromospheric 
dynamics. 3D simulations of solar-like stars including a chromosphere and 
corona are now computationally affordable (Leenaarts et al. 2007; Hansteen 
et al. 2007; Carlsson et al. 2010). 

On the observational side, a new generation of Fabry-Perot interferometers 
(FPI), for example IBIS at the Dunn Solar Telescope (DST) and CRISP at 
the Swedish 1-m Solar Telescope (SST), have provided evidences of very fine 
structure in the chromosphere. 

1.2.1 The chromospheric landscape 

The definition of chromospheric fine structure has evolved as new discoveries 
were made. It is widely accepted that the chromosphere includes the char- 
acteristic grass-like topology usually seen in H a images (Rutten 2006), but 
it is not clear where the boundaries of the chromosphere are. Fig. 1 .5 shows 
three Ca II 8542 A images acquired with SST/CRISP. This strong spectral line 
shows photospheric granulation at the wings and chromospheric fibrilar fea- 
tures in the core. Some chromospheric features are: 

• Straws are bright features seen in the core of chromospheric lines (Rutten 
2007). They start in facular regions in the photosphere and are much 
brighter than their surroundings in the chromosphere, showing hedge 
shapes in filtergrams. Fig. 1.5 shows a close view of straws, . In the upper 
panel, photospheric faculae become brighter than the surroundings at 
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Figure 1.5: Images in Call 8542 A taken at the SST. Top: the mid-photosphere. Mid- 
dle: High-photosphere/low-chromosphere. Bottom: the chromosphere. 
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mid-photospheric heights. Fibrils seem to originate in these straws, as 
shown in the lower panel in Fig. 1.5. 

• Fibrils, also known as mottles, are elongated dark features that in 
H I 6563 A (hereafter H a ) form a grass-like canopy covering internetwork 
cells at any heliocentric angle. They are also present Ca II images, but they 
only appear around network patches, as shown in Fig. 1.5. They are very 
dynamic and show transversal motions over a time period of 2 seconds. 
Hansteen et al. (2006) and Rouppe van der Voort et al. (2007) proposed 
that dynamic fibrils are driven by magneto- acoustic shocks that leak into 
the chromosphere along magnetic field lines. 

• Spicules: limb images of the chromosphere are dominated by spicules. De 
Pontieu et al. (2007) provided a detailed description of spicules and pro- 
posed physical mechanisms that could drive them. Spicules appear as thin, 
long highly dynamic features, usually reaching heights around 5000 km 
(see Fig. 1.4). Their width varies from 700 km down to current telescopes 
diffraction limit (~ 100 km). Spicules are classified in type I and type 
II. Type I spicules move up and down in time scales of 3-7 minutes and 
some of them present transversal motions. However, those fibrils that do 
not move transversely show acceleration and trajectories that are similar to 
those of dynamic fibrils, suggesting that they are also driven by magneto- 
acoustic shocks. Type II spicules are very dynamic and show apparent 
speeds between 50-150 km s , disappearing in time scales of 5-20 s. The 
mechanism driving type II spicules is not well understood, although their 
rapid disappearance suggests that strong heating could be ionizing Call 
atoms. Recently, Rouppe van der Voort et al. (2009) found on-disk coun- 
terparts of spicules, which produce a clear signature in the blue wing of the 
H I 6563 and Ca II 8542 lines. 

• Filaments are (dark) cold clouds of material that according to their 
temperature, belong to the chromosphere. They present typical lengths 
of 200000 km, with thickness's of 5000 km (Stix 2002). Towards the 
limb, filaments are known as prominences that appear hanging above 
the chromosphere up to 50 000 km. The only known mechanism that 
can sustain such cold and dense material is an electromagnetic force. 
Photospheric observations show that filaments predominantly exist along 
neutral magnetic field lines. Present observational efforts aim at measuring 
magnetic field in the filaments. 

1.2.2 Chromospheric heating 

One outstanding question about the chromosphere relates to its energy bud- 
get. Why is the outer atmosphere of the sun hotter than the photospheric sur- 
face? Semi-empirical, one-dimensional models of quiet Sun require the aver- 
age temperature to rise above the photosphere to reproduce the chromospheric 
intrinsic enhanced emission (VAL3, Vernazza et al. 1981). However, observa- 
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supergranulation Wedemeyer-Bohm el a!. BOW) 



Figure 1.6: Schematic structure of the quiet sun atmosphere (Wedemeyer-Bohm et al. 
2009), including some network patches. The black solid lines represent magnetic field 
lines that are anchored in the photosphere through network regions. 

tions show that the chromosphere is vigorously active and strongly inhomoge- 
neous. Carlsson & Stein (1995) demonstrated that enhanced emission can be 
produced by shocks without increasing the mean gas temperature. At present, 
the connection of waves with chromospheric heating appears widely accepted 
(Fossum & Carlsson 2005; Cauzzi et al. 2007; Wedemeyer-Bohm et al. 2007; 
Vecchio et al. 2009), however the exact role of these waves is still under de- 
bate. 

1.2.3 Magnetic field configuration 

The relatively organized and elongated fibrils seen in chromospheric lines (see 
Fig. 1.5) suggests the presence of magnetic fields. However, if magnetic fields 
dominate the chromosphere with j8 <C 1, those should be almost force free, 
leading to relatively smooth magnetic configuration. Thus, the complex fine 
structure must relate to the thermodynamics of the plasma (Judge 2006). 

The features described in §1.2.1 and their connection with the photosphere 
have been contextualized by Wedemeyer-Bohm et al. (2009) in the cartoon 
shown in Fig. 1.6. Magnetic field lines form a canopy where j8 ~ 1. Below 
the magnetic canopy, acoustic waves originating in the photosphere are 
dissipated, producing short-lived bright features seen in the wings of Call 
images. In this cartoon, fibrils and spicules form the magnetic canopy, 
which originates from network patches in the photosphere and extends over 
internetwork regions in the chromosphere. 



We investigate the relation between fibrils and magnetic fields in Paper III, 
looking for an observational evidence of their alignment. 
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2 Velocity references on solar 
observations 



This chapter describes the inherent difficulties in measuring absolute line-of- 
sight (LOS) velocities from spectroscopic observations. To illustrate the im- 
portance of this problem, we recall the discussion in Chapter 1 where we 
summarized the main theoretical frameworks that have been proposed to ex- 
plain the structure and dynamics of sunspots. A key difference between flux- 
tube models and the field-free gap model is that the latter explains penum- 
bral filaments as convective intrusions where magnetic field is weak enough 
not to suppress convection. Thus, observational evidence of convective mo- 
tions inside penumbral filaments would be a key to explaining the origin of 
its filamentary structure and choosing among existing models. These types of 
measurements are very hard to make because in the upper part of the filament 
overshooting convection is expected to be weak. At the same time, the pres- 
ence of the Evershed flow and the small scales involved, makes it very hard to 
establish the existence of overturning convection inside penumbral filaments. 
For these reasons, a very accurate velocity calibration is needed to make it 
possible to distinguish between the upflowing and downflowing components 
of any convection. 



2.1 The calibration problem 

The fundamental question that is addressed here is: What defines the local 
frame of rest on the Sun? An observer placed on the Sun would apply Eq. 2.1 
to measure LOS velocities, would apply the relationship 

v = ^ • c, (2.1) 

where A is the observed wavelength, A) is the reference wavelength (usually 
the laboratory wavelength of the line of interest) and c is the speed of light. 
However, ground-based observations are affected by the rotation of the Earth 
( v e,rot)> tne radial component of the Earth's orbital motion (v ffi or bit)> tne ro- 
tation of the Sun (vq ro t) and the graviational redshift (vgrav), the relarion is 
in reality more complex, 

A Ao 

v = — ^— • c + v ©,rot + v e:O rbit + v 0,rot + vgrav, (2.2) 
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Furthermore the precision of the atomic data limits the accuracy of the conver- 
sion from wavelength to velocities, regardless of how accurate the instrument 
is. This calibration issue becomes more severe when observational practical- 
ities make it difficult to compensate for the last three terms of Eq. 2.2. The 
gravitational redshift has a theoretical constant value of 633 m s _1 at the sur- 
face of the Earth (Cacciani et al. 2006). 

In addition, current instruments for solar observations seldom use labora- 
tory light sources as reference wavelengths. The obvious solution of finding 
a Ao on the Sun itself is confounded by the convective lineshifts discussed in 
Sect. 1.1.1. The magnitude of these shifts is different for different lines. 

Below, we summarize some of the methods that have been used to define a 
velocity reference for spectroscopic observations. 

1 . Convective motions in sunspot umbrae are suppressed by the presence of 
strong magnetic fields. Thus, it is common to assume that the umbra is at 
rest, defining a reference for line-of-sight velocities (e.g., Beckers 1977; 
Scharmer et al. 2008; Ortiz et al. 2010). However, this assumption usually 
does not hold higher up in the chromosphere, where umbral flashes associ- 
ated with shocks (Socas-Navarro et al. 2000a) produce strong blue-shifts. 
This approach carries the risk of being affected by spurious line shifts pro- 
duced by molecular blends that only form in the umbra because it is colder 
than the surrounding granulation. Eq. 2. 1 can be used to compute the con- 
version from the wavelength scale to a velocity scale, but in this case Ao 
is the central wavelength of the spectral line measured in the umbra of a 
sunspot. 

2. Telluric lines are sometimes present in the spectral range that is being ob- 
served. These lines are formed in the Earth's atmosphere. Thus, they al- 
low the definition of a very accurate laboratory frame of rest that can be 
converted to the solar frame using ephemeris constants, the time of the ob- 
servations and solar rotation (Eq. 2.2). The conversion from wavelength to 
velocities is calculated using the laboratory wavelength of the line of in- 
terest. Martinez Pillet et al. (1997) and Bellot Rubio et al. (2008) used this 
approach to calibrate their observations. 

3. A spectral atlas can be used to calibrate observations, as the effects of the 
rotation and translation of the Earth usually have been compensated for. 
Langangen et al. (2007) used the atlas acquired with the Fourier Trans- 
form Spectrometer at the McMath-Pierce Telescope (hereafter FTS atlas) 
of Brault & Neckel (1987) to calibrate some of his observations. This at- 
las was acquired at solar center, thus its usability is limited to disk center 
observations. 

4. Numerical models can be used to compute the convective shift of a line 
and use it as a reference for velocities. The advantage of this approach 
is that the convective shift is measured relative to the assumed laboratory 
wavelength of the line, so it is insensitive to uncertainties in the atomic 
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Table 2.1: Spectral lines used to create calibration data in Paper I. 



C I 


5380.34 


forms m deep layers of the photosphere 


re I 


5250.21 


magnetometry 


re I 


jZdv.oj 


magnetometry 


re I 




Doppler measurements 


re I 


6082.71 


abundance indicator 


re I 


6301.50 


magnetometry 


re I 


6302. 4y 


magnetometry 


re I 


/oyo.3o 


Doppler measurements 


Call 


8498.01 


chromospheric diagnostics 


Call 


8542.05 


chromospheric diagnostics 


Can 


8662.16 


chromospheric diagnostics 



data. Borrero & Bellot Rubio (2002) computed a two-components model 
from the inversion of photospheric Fel lines. This model only contains 
the vertical component of the velocity field, thus it is limited to solar cen- 
ter. Tritschler et al. (2004) and Franz & Schlichenmaier (2009) used this 
model to calibrate disk center observations. However, Bellot Rubio et al. 
(2004), used it together with the empirical results of Balthasar (1988) to 
estimate lineshifts also off solar centre. Langangen et al. (2007) used a 3D 
hydrodynamical simulation of solar convection to calibrate observations on 
C I 5380 A. The central wavelength of this line is not known with enough 
precision to allow the use of the FTS atlas method that was used with their 
observations. 

2.2 Calibration data from hydrodynamic granulation 
models 

In Paper I, we extend the calibration method employed by Langangen et al. 
(2007) for the C I 5380 line. Snapshots from a 3D hydrodynamical simulation 
are used to compute synthetic profiles assuming Local Thermodynamic Equi- 
librium (LTE) (see Fig. 2.1). The convective shift of the spatially-averaged 
profile is measured from spectra computed with the numerical simulation. 
Our calculations are performed for eleven selected lines of interest for solar 
observers (listed in Table 2.1), over a range of heliocentric angles (distance 
from solar disc centre). The synthetic line profiles are provided in digital form 
to the community. This method assumes that 3D models can reproduce the 
correlation between brightness and Doppler shift of granulation in a statistical 
sense (see Fig. 2.2). The elemental abundances is used as a free parameters 
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Figure 2.1: Left: Continuum intensity at 6301 A image from one of the 3D snap- 
shots from the numerical simulation. The latter has been resampled at lower spatial 
resolution. Right: Spectra from a blueshifted granule (solid black line) and redshifted 
intergranular lane (dashed line). The grey line corresponds to the spatial average of this 
snapshot. Two crosses on the left panel indicate the location of the spectra represented 
on the right panel. 
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Figure 2.2: Synthetic profiles and the corresponding bisectors resulting from the cal- 
culations of Fe I 5576.09 A. The grey scale indicates the variation of the heliocentric 
angle from jU = 1.0 to jU = 0.3. 



2.2 Calibration data from hydrodynamic granulation models 
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to achieve the best possible agreement between the synthetic profiles and the 
FTS atlas. The estimated parameters should not be regarded as abundances, as 
they also compensate for uncertainties in the atomic data, the LTE approxima- 
tion used in the radiative transfer calculations, and other errors. The accuracy 
of our results is inferred from experiments carried out using the 3D models. 
From this and from observational tests, we estimate the results to have an 
accuracy of 50 ms" 1 at solar disk centre. 

Our results have been analyzed using bisectors. The bisector of a spectral 
line indicates the center of the profile as a function of intensity. Fig. 2.2 illus- 
trate our results for the Fe I 5576.09 A line. The line bisectors are shown along 
with the profiles at different heliocentric angles. 

In addition to providing calibration data, Paper I discusses the variations 
of the bisectors with disk position (/i). This limb effect is found to mainly be 
caused by the 3D structure of the granulation, while the changing intensity- 
velocity correlation with height plays a minor role. 
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3 Chromospheric diagnostics 



It was mentioned in Chapter 1 that chromospheric observations are more dif- 
ficult than those of the photosphere, especially when polarimetric measure- 
ments are involved. Spectral lines that are sensitive to the chromospheric range 
are usually very broad and only the core, where less photons are emitted, 
shows chromospheric features (Cauzzi et al. 2008), as illustrated in Fig. 3.1 
where the granulation present in the wings smoothly changes into a chromo- 
pheric landscape close to the core of the line. The lack of light, in combination 
with a broad profile and weaker magnetic fields than in the photosphere con- 
spire to reduce the amplitude of Stokes Q, U and V profiles. 

The obvious solution to this problem would be to increase the exposure time 
of observations. However, the chromosphere is vigorously dynamic and long 
integration times usually translate into image smearing. The evolution time 
scale in the chromosphere can be estimated using an estimate of the Alfven 
speed v A {B = 100 G,z = 1000 km) ~ 10 5 m s 1 in the chromosphere (see 
page 83, Priest 1982). In the case of the SST, the diffraction limit at 854.2 
nm is 0."18 which corresponds to 130 km on the surface of the Sun. Thus, 
the chromosphere cannot be assumed to be static for times longer than 1.3 
seconds (see van Noort & Rouppe van der Voort 2006). 

Furthermore, this time scale also limits the spectral coverage of FPI obser- 
vations, as only one wavelength can be observed at the time. It is normally 
assumed that the Sun does not change during a full scan of the line. If the the 
spectrum is sampled using a large number of frequency points, this assump- 
tion may not hold. Therefore, observing the chromosphere involves a trade-off 
between sensitivity, cadence and wavelength coverage. 



3.1 Detectability of magnetic fields in the chromosphere 

During the past decade, the lines of the Call infrared triplet have been ex- 
tensively used to diagnose the chromosphere (see Langangen et al. 2008; 
Leenaarts et al. 2009; Cauzzi et al. 2009, and references therein), sometimes 
including polarization. Observational papers have usually studied cases with 
relatively strong magnetic field (Socas-Navarro et al. 2000a; Pietarila et al. 
2007a; Judge et al. 2010; de la Cruz Rodriguez et al. 2010), whereas theo- 
retical approaches have been restricted to ID models (Pietarila et al. 2007b; 
Manso Sainz & Trujillo Bueno 2010). 
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Figure 3.1: Observations of Call 8542 A. The images show a transition from the 
photospheric wings of the line to the chromospheric core. The dataset was acquired 
by Luc Rouppe van der Voort (ITA-UiO) with SST/CRISP. 



3.1 Detectability of magnetic fields in the chromosphere 
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Ca II Grotrian Diagram 
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Figure 3.2: Call atom model used to compute the calcium lines in Paper I & II. 

Transitions between bound states are marked with solid lines, whereas transitions from 
a bound state to a free state are represented by dashed lines. The wavelengths of the 
transitions are given in A. 



In Paper II a snapshot from a realistic simulation of the chromosphere 
is used for the first time to compute synthetic full Stokes spectra. We use a 
simplified Call model atom that consists of 5 bound levels plus ionization 
continuum, which is illustrated in Fig. 3.2. The populations of the atom are 
computed in non-LTE evaluating the 3D radiation field as in Leenaarts et al. 
(2009). 

This study is partially motivated by the ongoing debate on requirements of 
the instrumentation needed to observe chromospheric polarization in the quiet 
Sun using the Ca II infrared triplet lines. We study the combined effect of spec- 
tral resolution and noise on our simulated observations of the chromosphere, 
considering the following: 

• All the polarization is due to the Zeeman effect. We neglect the Hanle ef- 
fect, which depolarizes or polarizes the light depending on the scattering 
geometry and changes the ratio between Q and U (Manso Sainz & Trujillo 
Bueno 2010). 
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Figure 3.3: Spatially-averaged spectrum from the 3D simulation computed without 
microturbulence (solid-black line) and with microturbulence (solid-grey line). The so- 
lar FTS atlas is plotted for comparison (dashed-line). 

• The cores of our synthetic profiles are unrealistically narrow, probably be- 
cause the model is missing small scale random motions. Conclusions based 
on these results would underestimate the effect of noise and overestimate 
the effects of instrumental degradation. Thus, we use microturbulence to 
broaden our profiles to the same width that is observed in spatially-resolved 
profiles. In Fig. 3.3 the spatially-averaged spectrum from the 3D simulation 
with and without microturbulence are compared with the FTS atlas (see 
Brault & Neckel 1987). 

• Instrumental degradation is described by a Gaussian point spread function 
that operates on the spectra. Additive random noise following a Gaussian 
distribution is introduced after the convolution with the instrumental pro- 
file. 

Full Stokes monochromatic images computed from the 3D simulations in the 
Ca II 8542 A line are shown in Fig. 3.4. The images have a lot of sharp features 
that are partially produced by Doppler shifts of the line. As the chromospheric 
core of the synthetic spectra is unrealistically narrow and strong, intensity 
variation from Doppler shifts are stronger than in reality. 

Our results suggest that current FPI instruments are not sufficiently sensi- 
tive to detect circular polarization in the quiet Sun chromosphere using the 
Call 8542 Aline. 



3.2 Non-LTE inversions from a 3D MHD simulation 

Inversion codes have been extensively used to infer physical quantities from 
spectrometric and spectropolarimetric data. Inversions involve least-squares 
fits of the parameters of an atmospheric model, in order to reproduce observed 
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Figure 3.4: Stokes images at —75 mA from the core of the line. The Stokes Q and U 
images are scaled between ±0.002 and Stokes V is scaled between ±0.02, relative to 
continuum intensity. 



profiles. In the photosphere, LTE conditions or a Milne-Eddington atmosphere 
(Ruiz Cobo & del Toro Iniesta 1992; Bellot Rubio & Borrero 2002) are often 
assumed to simplify and speed up the radiative transfer calculations. Interest- 
ing results have been achieved with Milne-Eddington inversions of the chro- 
mospheric He I 10830 A lines (Lagg et al. 2004; Asensio Ramos & Trujillo 
Bueno 2009; Kuckein et al. 2010), which seem to form in the upper chromo- 
sphere (Centeno et al. 2008). 

The work presented by Socas-Navarro et al. (2000b) demonstrated that non- 
LTE inversions of solar observations are possible. Along those lines, Pietarila 
et al. (2007a) carried out inversions of Ca II 8542 A data, to measure chro- 
mospheric quantities in quiet Sun. More recently, de la Cruz Rodriguez et al. 
(2010) used the same scheme to carry out inversions on very high resolution 
observations of a sunspot showing umbral flashes. 

However, it is hard to quantify how accurately inversions can provide chro- 
mospheric information, with the commonly used assumptions related to the 
radiative transfer calculations: 

• The populations of the levels of the atom are computed in non-LTE assum- 
ing plane -parallel geometry. 

• Optionally, the computation of populations can be accelerated by neglect- 
ing the effects of the velocity field in the outcoming intensity. Thus, fewer 
azimuthal angles can be used to evaluate the radiation field. Under those 
conditions, the line profile becomes symmetric and only one half needs to 
be computed. 

• The fitted model is assumed to be in hydrostatic equilibrium to impose 
consistency between temperature and density. 
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Figure 3.5: Illustrative examples of the fits (solid-line) to the simulated observations. 
The left column is an example of a good fit, whereas the column on the right illustrates 
a poorer fit. The top row corresponds to Stokes / and the bottom row to Stokes V. 



In Paper II, we use the synthetic observations described in Section 3.1, with- 
out microturbulence, to test the Non-LTE Inversion Code based on the Lorien 
Engine (NICOLE) (Socas-Navarro et al. 2010). The results of the inversion are 
compared with the quantities from the 3D simulation model. The inversions 
provide a good estimate of the chromospheric average line-of- sight velocity 
and magnetic field. 3D non-LTE effects could be affecting temperature, which 
presents less temperature contrast than the original model. Fig. 3.5 shows two 
examples of fitted profiles from different pixels. The left panel corresponds 
to a good fit of the line, whereas the right panel shows a poor fit to the ob- 
served profiles. These failures originate the inversion noise that is mentioned 
in Paper II. 



3.3 Magnetic fields in chromospheric fibrils 

It is widely assumed that fibrils outline chromospheric magnetic fields. Fibrils 
usually appear around facular regions in Ca II 8542 filtergrams (Rutten 2007), 
supporting the connection between fibrils and magnetic fields. The goal is to 
find direct observational evidence of the alignment between magnetic fields 
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Figure 3.6: Superpenumbral fibrils in the surrounding of two sunspots observed in 
Ca II 8542 A. The left panel shows a wideband image of the photosphere whereas the 
panel on the right corresponds to a narrowband image acquired at -161 mA from the 
core of the line. Images acquired with the SST. 



and fibrils. We use two full Stokes datasets acquired in different telescopes 
with instruments of different type, to measure the orientation of magnetic 
fields in superpenumbral fibrils. The first dataset was acquired with SPINOR 
(Socas-Navarro et al. 2006) at the Dunn Solar Telescope (DST), a slit based 
instrument that allows a large wavelength coverage at a spatial resolution of 
0."6. The second dataset is acquired with SST/CRISP at very high cadence 
achieving a spatial resolution of 0."18 but with a limited spectral coverage 
(see Fig. 3.6). In these datasets, the Stokes Q and U spectra are integrated 
along the length of the fibrils in order to improve the S/N ratio. The azimuthal 
direction of the magnetic field (%) is calculated using the ratio between Stokes 
Q and U. 

tan(2 Z ) = ^ (3.1) 

Our measurements suggest that fibrils are mostly oriented along the magnetic 
field direction, however we find evidence of misalignment in some cases. This 
is both surprising, interesting, and hard to explain. Judge (2006) proposed that 
if j8 < 1 in the chromosphere, then magnetic fields should be almost force-free 
showing smooth spatial variations. The fine structure seen in chromospheric 
observations should then primarily be produced by the thermodynamic prop- 
erties of the gas. Our results could be compatible with this scenario. 
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4 Data collection and processing 



Some of the techniques that are described in this chapter are partially covered 
in Paper IV. However, the backscatter problem (see §4.3) and the telescope 
polarization model (§4.4) have not yet been described in separate publications, 
but are planned to appear in a forthcoming paper (de la Cruz Rodriguez et al. 
2011). 



4.1 The SST and CRISP 

The data presented in Paper III and Paper IV was acquired with the Swedish 
1-m Solar Telescope (SST) (Scharmer et al. 2003), located on the island of La 
Palma. Our narrow band data are acquired with the CRisp Imaging Spectropo- 
larimeter (CRISR Scharmer 2006) which is based on a Fabry-Perot interfer- 
ometer that allow for narrow band observations at very high spatial resolution 
and cadence, providing spectral information at the same time. Atmospheric 
turbulence is compensated for with adaptive optics (AO), in order to improve 
image quality. CRISP is mounted in the red beam of the SST (see Fig. 4.1). 
The light that has been corrected by the AO system passes through the chop- 
per and the pre-filter. Part of the light is reflected to the wideband camera. The 
other part is modulated with liquid crystals, producing linear combinations 
of the four stokes parameters. Afterwards, the light beam passes through the 
CRISP. The p and s polarizations are separated by a beam splitter into two 
beams that are recorded with separate cameras. 



4.2 Flat-fielding the data 

Science data taken with a CCD camera can be corrected for pixel-to-pixel in- 
homogeneities in the response of the camera. Normally, if the CCD is illumi- 
nated with a flat and homogeneous light source, intensity variations are mostly 
produced by pixel-to-pixel sensitivity variations, dirt and fringes. Thus, flat- 
field calibration images can be acquired to correct for these intensity varia- 
tions. 

A particular problem arises by the presence of time-dependent telescope 
polarization in the data. Normally, the flat-field images are taken at a different 
time than the science data. Thus, the amount of polarization introduced by 
the telescope can differ significantly between science and flat field data. In 
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Figure 4.1: Sketch of CRISP at the SST. The incoming light passes the chopper and 
the pre-filter that isolates the observed spectral line. A beam splitter separates 10 % of 
the light to the WB camera. The rest of the light passes through the liquid crystals and 
the Fabry-Perot etalons (HRE and LRE). Finally a polarizing beam splitter separates 
the p and s components of the polarized light to different cameras (1M). 



principle this should not be a problem if our cameras could detect Stokes 
parameters directly. Instead, four linear combinations of Stokes /, Q, U and V 
are acquired. If seeing were not present, the demodulation of the data could be 
carried out directly. However, in our case image restoration needs to be done 
in order to remove residual effects of seeing, not fully compensated for by the 
AO. As the image reconstruction is done with modulated data, artifacts appear 
if the flats and science data are not taken at similar times. 

Additionally, flat-fielding narrow-band images is more complicated than 
flat-fielding wideband images. In the case of CRISP, inhomogeneities on the 
surface of the FPI etalons produce field-dependent wavelength shifts of the 
transmission profile of the instrument. These are called cavity errors because 
they introduce variations in the FPI cavities that define the wavelength. The 
combined effect of cavity errors and the presence of a spectral line, produce 
field-dependent intensity variations purely introduced by the slope of the spec- 
tral line. At the same time, variations in the reflectivity of the etalons across 
the field-of-view translate into minor variations of the width of the instrumen- 
tal profile, and therefore also to overall transmission variations. This effect is 
much smaller than intensity fluctuations produced by cavity errors. Fig. 4.2 
illustrates these two effects on the Fe I 6302 A line. 

In Paper IV we address the flat-fielding problems produced by time varying 
telescope polarization and by cavity/reflectivity errors. We propose a method 
to flat-field polarimetric data affected by telescope polarization. A numerical 
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Figure 4.2: Left: A theoretical profile of CRISP is plotted with a grey line at two spec- 
tral locations. This example illustrates how a shift of the instrumental profile leads to a 
variation of the intensity on the flanks of a spectral line. The profile has been scaled to 
the intensity of the line at the location where it is centered. Right: Intensity variations 
produced by reflectivity errors. The line has been convolved with the instrumental 
profile of the same color. 



scheme is used to model and remove the fingerprints of the spectral line from 
our flat-field data. 



4.3 The backscatter problem 

The content of this section only applies to observations carried out at infrared 
wavelengths, and it was motivated by our first observations in Ca II 8542 A. 
The methods described here were developed together with Michiel van Noort 
and it is the main result of my first campaign with CRISP in 2008. 

The problem 

The CRISP acquisition system includes three back-illuminated CCD cameras 
(Sarnoff) that can record 35 frames per second. The quantum efficiency of 
these cameras decreases towards long wavelengths. Above 700 nm the CCD 
becomes semi-transparent, letting part of the light to pass through. Further- 
more, the images show a circuit-like pattern that cannot be removed by tradi- 
tional dark-field and flat-field corrections, as illustrated in Fig. 4.3. 
Examination of pinhole data and flat-fielded data shows that: 

1 . There is a diffuse additive stray-light contribution in the whole image. 

2. The stray-light contribution is much smaller in the circuit pattern and the 
gain appears to be enhanced. 
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Figure 4.3: Left: Subfield of an image that has been traditionally fiat-fielded, show- 
ing the circuit-like pattern. Right: The same subfield corrected using the numerical 
scheme that is proposed in this section. This colormap is chosen in order to enhance 
the contrast of fibrils and the electronic circuit. 

The model 

These problems can be explained by a semi-transparent CCD with a diffusive 
medium behind it, combined with an electronic circuit located right behind 
the CCD that is partially reflecting and therefore also less transparent to the 
scattered light. Fig. 4.4 represents the simplified structure of the camera. Un- 
der normal conditions an image recorded with a CCD, l Q , can be described in 
terms of the dark-field D, the gain factor Gf and the real image I r : 

I = D + (G r I r ) (4.1) 

However, in the infrared we need a more complicated model: 

Io = D + f(l-f)[(G b GfI r )*P]G h +fG f I r (4.2) 

where / represents the overall fraction of light absorbed by the CCD, Gb is 
the gain for light illuminating the CCD from the back which should account 
for the electronic circuit pattern. In the following, we refer to Gb as backgain. 
P is a Point Spread Function (PSF) that describes the scattering properties of 
the dispersive screen. In the backscatter term we assume that (1 — f)GfI r Gb 
is transmitted to the diffusive medium where it is scattered. A fraction of the 
scattered light returns to the CCD, passing again through the circuit. The car- 
toon in Fig. 4.4 shows a schematic representation of the structure of the cam- 
era and the path followed by the light beam. 

The numerical approach 

In order to obtain the real intensity I r , the PSF P, the back-gain Gb and 
the front-gain Gf must be known. The transparency factor is assumed 
to be smooth because the properties of the dispersive screen seem to be 
homogeneous across the field-of-view. This allows us to include / in the 
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CCD Circuit Dispersive screen 

Figure 4.4: Cartoon showing the conceptual model of the camera. The dark arrows 
indicate the incoming light beam from the telescope and the gray arrows represent the 
back-scattered light that returns to the CCD. 



front and the back gain factors 

Gy = y/l - f G b 
G f = fG f , 

so Eq. 4.2 becomes: 

I = D + Gy [(Gf/Gflr) * P] + Gfl r (4.3) 

This problem is linear and invertible, but a direct inversion would be expen- 
sive, given the dimensions of the problem. A numerical approach can be used 
to iteratively solve the problem. We define 

J=(G y G f ,I r )*P (4.4) 

In the first iteration, we initialize / assuming that the smearing caused by P 
is so large that / can be approximated by the average value of the observed 
image multiplied with the back gain. The back gain Gy is assumed to be 1 for 
every pixel in the first iteration. Furthermore, we assume that the product Gfl r 
can be estimated from Eq. 4.1 ignoring back scattering, i.e., Gfl r I Q — D. 
These values are of the same order of magnitude as the final solution, and 
therefore correspond to a reasonable choice of initialization. 

Jo ~ (GyG f I r ) « (Gy (I —D)), (4.5) 

where the initial guess of the PSF P represents an angular average obtained 
from a pinhole image. The small diameter of the pinhole only allows to esti- 
mate accurately the central part of the PSF. The wings of our initial guess are 
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extrapolated using a power-law. The estimate of Gfl r is then 

G f I r = I -D-GyJ (4.6) 

which can be used to compute a new estimate J. This procedure is iterated, 
until G fl r and I are consistent. This new value of / is used to improve our 
estimate of Gy , by applying the same procedure to images that contain parts 
being physically masked (I r = 0), as in Fig. 4.5. In the masked parts where 
I r = we have, 

I o -D-GyJ = 0, (4.7) 
so the back-gain can be computed directly: 

Gy = 1 ^^. (4.8) 

Thus, every pixel must have been covered by the mask at least once in a cali- 
bration image in order to allow the calculation of the back-gain. With the new 
Gy we can recompute a new estimate of Gfl r . 

We now need to specify a measure that describes how well the data are 
fitted by the estimate of P and Gy(P). Since both Gfl r and Gy are computed 
based on self-consistency, we iteratively need to fit only the parameters of the 
PSF. We assume that the PSF is circular-symmetric and apply corrections to 
the PSF at node points placed along the radius. 

We use Brent's Method described by Press et al. (2002) to minimize our 
fitness function. This algorithm does not require the computation of deriva- 
tives with respect to the free parameters of the problem. When the opaque 
bars block a region of the CCD, an estimate of the back gain can be calculated 
for a given PSF according to Eq. 4.8. In our calibration data, the four bars of 
width L are displaced 0.5 L from one image to the next (see Fig. 4.5). This 
overlapping provides two different measurements of the back gain on each re- 
gion of the CCD. However, as the location of the bars changes on each image, 
the scattered light contribution is different for each of these measurements of 
the back gain. Our fitness function minimizes the difference between these 
two measurements of the back gain. 

Having thus obtained the PSF P and the back-gain Gy, we obtain Gf by 
recording conventional flats and assuming I r is a constant in order to obtain 
Gf fromEq. 4.9. 

Results 

The numerical scheme described in §4.3 produces the backgain (Gy) and an 
approximate PSF (P) that describes the scattering problem. The results are 
illustrated in Fig. 4.6. The back gain image shows a background with vertical 
dark areas in the extended gaps where the circuit pattern is not present. Those 
hollows probably indicate that the PSF is not perfect, an expected result given 
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Figure 4.5: Calibration images used to infer the scattering PSF and the backgain of 
the camera. The bars are displaced image-to-image so every pixel is covered at least 
once by the black bars. The images are shown with logarithmic intensity scale. 




Figure 4.6: Left: Fitted back-gain for one of the cameras. Right: Circular-symmetric 
PSF obtained from our fitting routine. The PSF is displayed using a logarithmic scal- 
ing. 
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Figure 4.7: Left: Values of the back gain along the x— axis at Y = 450. The black 
markers indicate the location of the hollows resulting from the limited accuracy of the 
PSF. Right: The radial variation of P, represented in logarithmic scale. 



the assumptions imposed on the shape of the PSF: it is constant over the entire 
field of view, and only radial variations are allowed. The PSF has extended 
wings that complicate the convergence of the solution, given the nature of 
our calibration data: the spacing between bars in the horizontal direction is a 
limiting factor to constrain such extended wings. 

The flat-fielding method described in the previous sections has been imple- 
mented in the image reconstruction code MOMFBD (van Noort et al. 2005) 
and is only used when a back-gain and a PSF are provided. Eq. 4.9 is reordered 
to correct the science data. 

Gfl r =I„~ Gy \{GyG f I r ) * P] - D (4.9) 

Since G fl r appears in both sides of the equation, some iterations are needed 
to estimate the term Gy[(GyG fl r ) *P], which is a computationally expensive 
process if thousands of images are flat-fielded. 

In Fig. 4.3 we show a frame that has been traditionally flat-fielded (left 
panel) and the same image flat-fielded according to Eq. 4.9. The method de- 
scribed in the present work removes the electronic circuit pattern from the 
images. The assumption imposed on the PSF allows us to estimate the back 
gain with limited accuracy, as is obvious from the Fig. 4.7. The uncertainties 
present in the PSF and the back gain are likely to affect the contrast of the 
corrected images. 
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Figure 4.8: The 1-m polarizer mounted on the entrance lens of the SST. Calibration 
datasets were acquired during entire days to characterize the polarimetric properties 
of the telescope. 



4.4 Telescope polarization model at 854.2 nm 

The turret of the SST contains optical elements that polarize the incoming 
light. Selbing (2005) studied the polarizing properties of the telescope at 630.2 
nm and proposed a theoretical model to characterize its temporal variation. 
Calibration images have been taken using a 1-m polarizer mounted at the en- 
trance lens of the SST (see Fig. 4.8). The polarizer rotates 360° in steps of 5° 
and several frames are acquired on each polarizer angle. Data were adquired 
during the whole day. These data were used to determine the parameters of 
the model proposed by Selbing (2005) at 854.2 nm. 

Telescope model 

Each polarizing optical element of the telescope is represented with a Mueller 
matrix. Mueller matrices of mirrors are noted with M and have two free pa- 
rameters. Assuming a wave that propagates along the z axis and oscillates in 
the x — y plane, the parameters are the de- attenuation term (R) between the 
xy components of the electromagnetic wave and the phase retardance (8) pro- 
duced by the mirror (see Selbing 2005). This form of M assumes that Q is 
perpendicular to the plane of incidence. In the more general case, M can be 
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written as, 
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R(a) corresponds to a rotation to a new coordinate frame, rotated an angle a. 
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L represents the entrance Lens. The form of the lens matrix represents a com- 
posite of random retarders, so it modulates the light without (de)polarizing it. 
The reference for Q is aligned with the lm linear polarizer axis and the values 
of the matrix are measured in that frame. 
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The model is built using the Mueller matrix of each polarizing element in the 
telescope (Eq. 4.10), as a function of the azimuth (<p) and elevation (6) angles 
of the Sun at the time of the observation. We have included the conversion 
factor from degrees to radians in the matrices, thus all the angles are given in 
degrees. 
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where, 

• co, ci, C2, C3 and C4 are the parameters of the entrance lens. 

• C5 and C(, are the de-attenuation and phase difference of the azimuth mirror. 

• C7 and eg are the de- attenuation and phase difference of the elevation mirror. 

• eg and cio are the de-attenuation and phase difference of the Schupmann 
mirror. 

• c\\ is the angle of the field mirror. 

We used a Levenberg-Marquardt algorithm (see Press et al. 2002) to fit the 
12 parameters (c) of the model to the calibration data. The orthogonality of 
1-m polarizer states is maximum every 45°. Thus only four angles are used in 
our fitting routine: 0°,45°,90°, 135°. The quality of the fit does not improve 
substantially by including data from other angles. The derived parameters of 
the Mueller matrix of the telescope are shown in Table 4. 1 . Unfortunately, we 
have not been able to analyze the errors due to time constraints. 

The quality of the calibration data is limited by the quality of the 1-m linear 
polarizer. The problem is posed in such a way, that we cannot measure the 
extinction ratio of the polarizer and the parameters of the lens at the same time. 
In our case, the polarizer has a significant leak of unpolarized light at 854.2 
nm. Using small samples of the sheets used to construct the 1-m polarizer, we 
have estimated the extinction ratio of the 1-m polarizer to be approximately 
0.4, so the parameters of the lens have been determined assuming that value. 
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Data collection and processing 



Table 4.1: Parameters of the telescope model obtained from calibration data. 
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In Fig. 4.9 the time dependence of the Mueller matrix of the telescope is 
shown along a whole day. The largest changes occur when the sun is close to 
zenith, because of the rapid movement of the telescope. The linear polarization 
reference is defined by the first mirror after the lens, however this is not very 
useful in practice because the turret introduces image rotation along the day. 
Instead, we use the solar north as a reference for positive Q by applying an 
extra rotation to M (e /. The rotation is produced by reflections inside the turret 
and by the variation of the angle between the first mirror after the entrance 
lens and the solar north along the day The angle between the first mirror and 
the solar north (j8) is computed by the telescope software every 30 seconds. 
Eq. 4.11 transforms the reference of Q and U to solar North-South axis. 

M' tel =M tel -R(P) (4.11) 

The only remaining question is the location of the solar north in our science 
images. At this point, Stokes Q and U are relative to the solar North-South 
axis, but there is no coupling between the polarization calibration and the 
image orientation. The angle between solar north and the horizontal on the 
optical table is 

ao = <p-e-TC-p, (4.12) 

where (p is the azimuth, 6 is the elevation, TC is the table constant and /3 is 
the tilt angle between first mirror in the telescope and solar north. The table 
constant is relative to the orientation of the optical table and it is +48° for the 
current setup. 

Fig. 4.10 and 4.11 show monochromatic Stokes /, Q, U, V images ac- 
quired in Ca II 8542 A. The dataset has been restored using the de-scattering 
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Figure 4.9: Mueller matrix of the telescope computed for 2009-05-28. Each panel 
corresponds to one element of the matrix placed in the same order. The solid line 
represents the time dependence of each value. The dashed line is the same value mul- 
tiplied by a factor 10. 





Figure 4.10: Stokes / (top) and Q (bottom) monochromatic images acquired in Call 
8542 A at -161 mA from the core of the line. Stokes Q is scaled to ±3% of the 
continuum intensity. 
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Figure 4.11: Stokes U (top) and V (bottom) monochromatic images acquired in Ca II 
8542 A at —161 mA from the core of the line. The images are scaled to ±3% of the 
continuum intensity. 
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Data collection and processing 



scheme from §4.3 and the telescope model presented in the current section. 
This dataset is used in Paper III to measure the alignment between fibrils and 
magnetic field. 
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5 Summary of papers 



Paper I: Solar velocity references from 3D HD 
photospheric models 

The aim of this paper is to help observers to accurately calibrate line-of-sight 
velocities. We use a 3D hydrodynamic simulation of the solar photosphere to 
compute spatially-averaged spectra that can be used as absolute velocity ref- 
erences. The line profiles are computed at different heliocentric angles, from 
disk center towards the limb. 

Our synthetic profiles are compared with observational data and several ex- 
periments are computed to estimate the accuracy of our method, which has an 
estimated error of approximately ±50 m s _1 at disk center. Our tests suggest 
that the variation of the bisectors towards the limb, is mostly produced by the 
3D topology of the photosphere. 

In Paper I, I carried out all the calculations and prepared all figures. The 
collaborators contributed to the scientific discussion and assisted in the writ- 
ing. 

Paper II: Non-LTE inversions from a 3D MHD 
chromospheric model 

In Paper II, we create synthetic full-Stokes observations in Ca II 8542 A from 
a snapshot of a realistic 3D simulation of the solar atmosphere. These observa- 
tions are used to estimate the amplitude of the Stokes profiles in quiet Sun. We 
discuss the effect that spectral degradation and noise have on our observations 
and discuss possible requirements for future instrumentation. 

In the second part of the paper, we use our synthetic observations to test 
our non-LTE inversion code. The fitted model is compared with the quantities 
from the 3D snapshot. We conclude that the inversion code is able to estimate 
the average chromospheric value of magnetic field, line-of-sight and veloc- 
ity. 3D non-LTE effects seem to affect the fitted temperature that in general 
presents less contrast than the original model. 

My contribution to Paper II was to compute the full-Stokes simulated ob- 
servations using the population densities provided by J. Leenaarts and carried 
out the Non-LTE inversions of the data. For that purpose I wrote an improved 
parallel version of NICOLE using MPI and a master-slave scheme. I prepared 
the main structure of the text in the paper and created all the figures. 
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Summary of papers 



Paper III: Are solar chromospheric fibrils tracing the 
magnetic field? 

The aim of this letter is to obtain an observational evidence that confirms 
the alignment between chromospheric fibrils and magnetic field. We use two 
datasets acquired with SST/CRISP and DST/SPINOR to measure the orienta- 
tion of magnetic field along chromospheric fibrils. We find that many fibrils 
are aligned with magnetic field, however in both datasets there are evidences 
of misalignment in some cases. 

For this paper I provide a restored dataset from CRISP, compensated for 
telescope polarization and with a calibrated reference for Stokes Q and U. My 
co-author prepared the SPINOR data and contributed to the scientific discus- 
sion and the writing. 

Paper IV: Stokes imaging polarimetry using image 
restoration at the Swedish 1-m Solar Telescope II: A 
calibration strategy for Fabry-Perot based instruments 

The image restoration step that is applied to our data, decouples the 1-to-l 
relation between the pixels of the CCD. When image reconstruction is applied 
to data showing sharp spatial variations produced by intrumentation, artifacts 
can appear. We propose a flat-fielding scheme for polarimetric data acquired 
with CRISP. We discuss the effect of the polarization introduced by the tele- 
scope and the optical setup in our flat-field data. In order to correct for spurious 
intensity fluctuations from cavity errors and reflectivity errors, we use a nu- 
merical framework that allows to model the spectral line on each pixel of the 
CCD. 

My contribution to this paper was to implement the numerical scheme that 
is used to model the flat fields and remove the intensity fluctuations produced 
by cavity errors and reflectivity errors. I contributed to the scientific discussion 
and prepared some of the figures in the paper. 
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6 Publications not included in this 
thesis 



• CRISP Spectropolarimetric Imaging of Penumbral Fine Structure 

Scharmer G. B., Narayan G., Hillberg T., de la Cruz Rodriguez J., Lof- 
dahl M. G., Kiselman D., Sutterlin P., van Noort M., Lagg A., 2008, ApJ, 
689, L69. 

• The magnetic SW Sextantis star RXJ1643.7+3402 

Rodriguez-Gil P., Martmez-Pais I. G., de la Cruz Rodriguez J., 2009, MN- 
RAS, 395, 973. 

• High-order aberration compensation with Multi-frame Blind Deconvo- 
lution and Phase Diversity image restoration techniques 

Scharmer G. B., Lofdahl M. G., van Werkhoven T. I. M., de la Cruz Ro- 
driguez J., 2010, A&A, 521, A68 

• Observation and analysis of chromospheric magnetic fields 

de la Cruz Rodriguez J., Socas-Navarro H., van Noort M., Rouppe van 
der Voort L., to appear in Proceedings of the 25th NSO Workshop: 
Chromospheric Structure and Dynamics, Memorie della Societa' 
Astronomica Italiana. 
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